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We numerically study models of H 2 and LiH molecules, aligned collinearly with the linear polarization of the 
external field, to elucidate the possible role of correlation in the enhanced-ionization phenomena. Correlation is 
considered at different levels of approximation with the time-dependent generalized-active-space configuration- 
interaction method. The results of our studies show that enhanced ionization occurs in multielectron molecules, 
that correlation is important and they also demonstrate significant deviations between the results of the single- 
active-electron approximation and more accurate configuration-interaction methods. With the inclusion of corre¬ 
lation we show strong carrier-envelope-phase effects in the enhanced ionization of the asymmetric heteronuclear 
LiH-like molecule. The correlated calculation shows an intriguing feature of cross-over in enhanced ionization 
with two carrier-envelope-phases at critical inter-nuclear separation. 

PACS numbers: 31.15.-p, 32.80.Fb, 33.80.Eh, 42.50.Hz 


I. INTRODUCTION 


When a diatomic molecule is exposed to ultra-short laser 
pulses, the ionization yield is enhanced at certain critical inter- 
nuclear distances. This phenomenon is known as enhanced 
ionization (El) and has been observed both experimentally [1- 
7] and discussed theoretically [8-13] (see Ref. [14] for a re¬ 
view on El in molecules). The main interesting feature that 
distinguishes ionization of a diatomic molecule from ioniza¬ 
tion of an atom is the double-well structure of the potential. 
At the equilibrium inter-nuclear separation, the molecular po¬ 
tential is dominated by the Coulombic monopole and in this 
sense it is similar to the atomic potential. As the inter-nuclear 
separation R increases, the internal barrier between the two 
wells broadens. In the presence of a strong electric field the 
diatomic potential is distorted and the potential changes dras¬ 
tically. Since the barrier height in one of the well reduces, the 
electron can directly tunnel to the continuum and this effect 
enhances the ionization probability at certain critical inter- 
nuclear distance. For one-electron systems, the localization 
of the wave function in the upper well of the effective poten¬ 
tial formed by the molecular potential and the external field 
is a key mechanism as explained in Ref. [8], It was shown 
earlier that electron localization in the presence of a strong 
field is an often encountered phenomenon [15, 16]. The prob¬ 
ability of localization in one of the wells depends on the laser 
intensity and frequency as well as on the phase [16], It is 
also important to emphasize that the localization of the elec¬ 
tronic wave packet depends on the preparation of the initial 
state. In Ref. [9] a different mechanism for El was proposed. 
In that work El is seen as the results of the creation of a pair 
of charge-resonant (CR) states which strongly couples to the 
external field at the critical inter-nuclear distance. Along with 
the creation of CR states the potential experienced by the elec¬ 
tron is also changed in the presence of the strong field. This is 
also known as charge-resonance-enhanced-ionization (CREI) 
of diatomic molecules. A one-dimensional model of H 2 was 
studied in Ref. [17] and a full three-dimensional calculation 


of H 2 was presented recently [18], For a one-electron sys¬ 
tem like H,j~, it has been widely discussed from which well of 
the diatomic potential the electron will escape. In Ref. [11], 
it is argued that the ionization can occur by resonant transi¬ 
tion to an excited state supported by the lower well. Using 
H 2 , it is also argued that the El in diatomic molecules occurs 
due to level crossings between the ground state and that of the 
excited states which dissociate into ionic fragments [19, 20], 
The El in more complex, e.g., asymmetric molecules was 
studied in Ref. [21-23], Experimentally the El was found in 
molecular iodine [1], and in molecular ions [2, 3], The El 
in CO 2 was experimentally studied in Ref. [4], and in Ref. 
[5] the tunneling site of the electron that takes part in the El 
of molecules was experimentally probed. Recent experiments 
directly detected El in CO and N 2 [6] and demonstrated a two 
peak structure of El in 11 2 [7], 

Solving the time-dependent Schrodinger equation (TDSE) 
directly in the presence of a strong external field is possi¬ 
ble only for atoms and molecules with one or two electrons. 
Most of the previous calculations therefore use the single- 
active-electron (SAE) approximation. They may, therefore, 
neglect some important features which may arise due to 
electron correlation. The purpose of this work is to elu¬ 
cidate El further by investigating in a systematic manner 
the role of electron correlation on this mechanism. Sev¬ 
eral time-dependent many-body methods have been devel¬ 
oped to take into account electron correlation as accurately as 
possible. The time-dependent /-(-matrix approach [24-27], 
the time-dependent configuration-interaction singles (TD- 
CIS) approximation [28-31], the multiconfigurational time- 
dependent Hartree-Fock method [32-38] and its generaliza¬ 
tions [39-43] have been applied to processes involving laser- 
matter interactions. The time-dependent restricted-active- 
space configuration-interaction (TD-RASCI) method [44], 
and time-dependent generalized-active-space (TD-GASCI) 
configuration-interaction method [45, 46] systematically take 
into account the electron correlation through a configuration- 
interaction expansion within a mixed-basis set approach with 
localized and grid-based orbitals to allow for an accurate 
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description of both, the bound and continuum states. This 
method includes the limits of full Cl, the SAE and the CIS 
approximations. The TD-GASCI method has been used to 
calculate the photoelectron spectra of one-dimensional two 
and four electron atoms and molecules [45] and the TD- 
RASCI method to calculate photoionization cross sections of 
beryllium and neon [44], In the present work, we use the 
TD-GASCI method to study correlation phenomena in El of 
model molecules. First, we consider a model of H 2 . In this 
homonuclear system we study the convergence of the TD- 
GASCI method by comparing it with fully correlated (“ex¬ 
act”) TDSE calculations. To study the correlation trend we 
use linearly polarized laser fields polarized along the inter- 
nuclear axis within the fixed-nuclei approximation. Further 
we consider a model of LiH and report on El phenomena 
within four electron molecules. In this asymmetric system 
we study the correlation effects in El and in addition the sen¬ 
sitivity to the relative orientation between the molecule and 
the dominant half-cycle of the short pulse as controlled by the 
carrier-envelope-phase (CEP). 

The paper is organized as follows. In Sec. II, we briefly 
present the TD-GASCI method and the special basis set used. 
In Sec. Ill, we use the TD-GASCI method to calculate the 
ionization yield as a function of inter-nuclear distance to study 
the El phenomena. We demonstrate the accuracy of the TD- 
GAS-CI results of the H 2 model by comparison with TDSE 
results. Further we consider calculations for the LiH model 
with different GAS partitions corresponding to accounting for 
electron correlation at different levels of approximation. In 
Sec. IV, we summarize our work and conclude. 


II. THEORY 


Here we briefly recall the TD-GASCI methodology [45]. 
To ease notation we consider a diatomic molecule with nuclei 
at ±R/2 and charges Z\ and Z 2 * respectively. The TDSE 
for N e i -electrons with fixed nuclei then reads (we use atomic 
units throughout this work) 

= m\m) ■ (i) 

The time-dependent Hamiltonian is given by 

N el N Bl 

m = E hi (f ) T ^ ^ tn ij , (2) 

i=l i<j 

with the one-body part of electron i given by 


hi(t) 



Z 2 


+ Vi ■ F(t) , (3) 


with F(t) is the laser field. In Eq. (2) the two-body interaction 
is given by 


Wij = 


ri - r. 


( 4 ) 
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FIG. 1. Schematic of the partition-in-space concept. The molecule 
with binding length R is centered on the grid, the binding potential of 
a LiH-like model is indicated by v(x) (dashed line). The basis set is 
partitioned into a central part |a;| < x c and an outer part. For the for¬ 
mer, localized Hartree-Fock and pseudo-orbitals for the description 
of the bound electrons are constructed. Details are given in Ref. [45]. 


The many-body wave function is expanded into a basis of 
Slater determinants |$/), 

!*(*)> = E C M*i) . ( 5 ) 

/evE*c 

with time-dependent expansion coefficients C’j(i) = 
(<l>/|\l/(f)) and I is a multi-index specifying the configurations 
drawn from the Hilbert space Vexc- The Slater determinants 
themselves are constructed from Nt> time-independent single¬ 
particle orbitals <pi (2A), spin orbitals), within a mixed-basis 
set approach. Thereby, the space variable r is partitioned into 
a central region M < |t c |, where localized Hartree-Fock (HF) 
and pseudo-orbitals are used to describe the bound electrons, 
and an outer region |r c | < |r| < |r e | with a finite-element 
discrete-variable-representation (FE-DVR) basis description 
for the outgoing electron wave, see Fig. 1 for a schematic of 
the ID situation (x corresponds to one component of r). De¬ 
tails of the basis and its implementation are given in Ref. [45]. 

The corresponding matrix form of the TDSE then reads as 

l Wt Cl{t)= ^ (6) 

J £ Vexc 

with the Hamiltonian matrix Hjj(t) = ( < I>/|TT(i)| < f > j). This 
matrix can be constructed by evaluating and partially trans¬ 
forming the associated one- and two-electron integrals. If the 
sums in Eqs. (5) and (6) run over all possible excitations Vk xc , 
the full Cl (FCI) method [47] is obtained and an exact solution 
of the TDSE within the limits of the finite single-particle ba¬ 
sis set is retrieved. Typically, an FCI expansion of the TDSE 
is numerically intractable for processes involving one or more 
continua as in photoionization, and for systems consisting of 
more than two active electrons due to an exponential scaling in 
the number of configurations with the number of basis func¬ 
tions. The generalized-active-space (GAS) concept aims to 
overcome this exponential barrier by choosing the most rele¬ 
vant Slater determinants for the dynamics under consideration 
and thus forming a subset of the FCI many-particle basis set, 
Vexc = V’gas in Eq. (5). By systematically including more 
Slater determinants, convergence of the method towards the 
fully correlated result is obtained. 

The relevant GAS partitions for this work are summarized 
in Fig. 2. The simplest approximation is the SAE approxima¬ 
tion, in which only one selected electron may occupy excited 
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FIG. 2. (color online). Generalized active spaces used in this work. 
The energies of the single-particle levels (up to N b spatial orbitals) 
are denoted by Ei. The two spin configurations, | J |') and | J_), are 
degenerate in the spin restricted HF ansatz. The left panel shows the 
SAE approximation, where only one electron (]. in level E 2 in this 
case) is allowed to be excited. The right panel shows the complete- 
active space (CAS*) situation, where correlation can be systemati¬ 
cally built in with the incorporation of more and more active orbitals. 
N e i -electrons are allowed to occupy K spatial orbitals and one elec¬ 
tron excitations for any spin configuration are allowed above level 
Ek- 
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FIG. 3. (color online). Normalized electric field with (a) i/jcep = 0 
and (b) </j C ep = n. 


Slater determinants, likewise the CIS approximation can be 
obtained from the GAS. By construction, no additional poten¬ 
tials are needed to be evaluated to obtain this widely used trun¬ 
cation. The other partition refers to the complete-active-space 
(CAS) concept [48, 49], which corresponds to a FCI descrip¬ 
tion (with possibly a frozen core) up to a spatial orbital index 
K. In the present TD-GASCI approach, the CAS is accompa¬ 
nied by additional single excitations out of the CAS, indicated 
in our notation by CAS* (N e j, K), where N e i denotes the ac¬ 
tive electrons and I\ the number of spatial orbitals within the 
CAS. In this context we wish to emphasize that we consider 
only one electron in the continuum throughout and double ion¬ 
ization is neglected, per construction of the GAS [45]. 


III. CORRELATION EFFECTS IN ENHANCED 
IONIZATION 


To study EI in molecules, we consider colinear models of 
H 2 and LiH. We use the regularized Coulomb potential be¬ 
tween the electrons and the nuclei. 



Here x is the electron coordinate, R is the inter-nuclear dis¬ 
tance, and Zi (1 = 1, 2) the nuclear charges. For all calcula¬ 
tions we use the softening parameter, s = 1. A two-electron 
model of H 2 is defined by Z\ = Z 2 = 1 and a four electron 
model of LiH by Z\ = 1, and Z 2 = 3. The electron-electron 


interaction in this model reads 


w(Xi,Xj ) 


1 

■\/{.Xi - Xj ) 2 + s 2 


( 8 ) 


Such kind of models for diatomic molecules are well 
established in strong-field physics involving near-infrared 
fields [50]. In the calculations we use the length gauge to 
describe the interaction with the laser field. We use the vector 
potential with a sine-square envelope [51], 

A(t) = — sin 2 sin (cut + <A:ep) (0 < t < T) , 
uj V 1 / 

(9) 

where T = n^j- is the pulse duration with n = 1, the number 
of cycles. The corresponding electric field is F(t) = — . 

F[) is the maximum amplitude of the laser field to be speci¬ 
fied below in the calculations, w is the angular frequency and 
(pCHP is the carrier-envelope-phase (CEP). The form of the 
laser pulse is shown in Fig. 3. For the present work we use 
800 nm (u> = 0.057) and single-cycle pulses of duration 2.65 
fs (110 a.u.). The vector potential in Eq. (9) ensures that the 
time-integral over the electric field is vanishing, as it should 
be [52]. We note that the pulse duration considered here is 
so short that a molecule will not have time to dissociate to 
the range of critical inter-nuclear distances for EI during the 
presence of the pulse. Experimentally a pump-probe scheme 
is therefore required to investigate the EI phenomena for such 
short pulses. 

To extract the total ionization probability, we add a complex 
absorbing potential (CAP) to the Hamiltonian, 

H CAP (t) = H(t) - iV cap , 


( 10 ) 
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where the CAP has the form [53] 


VcapO) = 1 - cos 


/ 7r(|x| - x C ap) \ 
\ 2(x s - Xcap) ) 


( 11 ) 


for |x| > .x'cap and zero otherwise. For sufficiently long prop¬ 
agation time, tf, after the end of pulse, the continuum part of 
the wave function has passed into the region of space beyond 
xcap and has been absorbed. The total ionization yield is then 
given by 


= (i2) 

withA/"(f/) = [54]. To extract the ionization 

yield after the end of the pulse, we propagate the equations of 
motion to a final time, tf = 241 fs. We found that this time 
is sufficient to obtain converged results, also for correlated 
situations. 

We partition the full simulation box into an inner (central) 
and an outer region. In the present simulation the main focus 
is on a wide range of inter-nuclear separations. Therefore, we 
need a relatively large box for the central region, \x\ < x c , 
as the correlated bound part of the wave function is described 
by localized HF and pseudo-orbitals. We use x c = 16 in the 
present calculation, cf. Fig. 1, and checked carefully for con¬ 
vergence of our results. We use 16 elements which have 8 
DVR functions each for the central region which results in to¬ 
tal 111 FE-DVR functions in the central region for all calcula¬ 
tions. Further increasing the number of DVR functions leads 
to almost no change in the HF energy. The TD-GASCI calcu¬ 
lation is done on the full simulation box to take into account 
the continuum. In the present calculation we set the box size 
to x e = 200, and the outer region is described by the FE-DVR 
functions. By using the same FE-DVR parameters as for the 
inner region, we therefore include 699 FE-DVR functions for 
the whole simulation box. We checked the convergence with 
respect to different box parameters and the above parameters 
gave converged results. 


A. Two-electron fD-like molecule 
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FIG. 4. (color online). Total ground-state energy of the 2-electron 
H 2 -like molecule. In the upper curve the SAE and the CIS approx¬ 
imation overlap with the HF result. CAS* (2, 2) represents the TD- 
GASCI calculation with two active orbitals in the GAS. Similarly 
CAS* (2,11) represents the TD-GASCI calculation with eleven ac¬ 
tive orbitals in the GAS. 


coefficients. In this approximation the sum runs over all vir¬ 
tual orbitals, a, with a fixed core orbital, i. Therefore, it rep¬ 
resents the effective interaction felt by the electrons from the 
ionic core. In the CIS approximation the wave function reads 


|4'cis(i))=Co(f)|4>o)+^C'“(f)|$“) . (14) 

i,a 


The two electron system is the preferred choice of system to 
use to validate the approach because we can directly compare 
the results of the TD-GASCI simulations with TDSE simu¬ 
lations, which are performed in the same FE-DVR basis set. 
We begin with the ground state energy calculation with imag¬ 
inary time propagation to demonstrate the features of the TD- 
GASCI method. In Fig. 4 we show the total ground state en¬ 
ergy of ID H 2 with HF, SAE, CIS approximations and differ¬ 
ent GAS methods along with the TDSE calculations. Let us 
first discuss the approximations including singly excited de¬ 
terminants only. In the SAE approximation, the wave function 
reads 

I’Wf)} = Co(i)l^o) + E <?“(*)I*“> • (!3) 

a 

Here, |$ 0 ) i s the HF reference determinant, and <[>'') is 
a singly-excited determinant with the time-dependent CI- 


This approximation incorporates all possible single excita¬ 
tions from the core orbitals to the virtual orbitals. The time- 
dependent CIS approximation has been widely used for under¬ 
standing the processes involving strong-held ionization [29- 
31]. Although these two approximations take into account 
a certain contribution to the correlation by accounting for the 
singly-excited determinants, the SAE and CIS approximations 
do not improve the ground state energy. This is according 
to Brillouin’s theorem [47], which states that the matrix el¬ 
ements of the Hamiltonian between the reference state and 
singly-excited states vanish within a HF basis. From Fig. 4 
we can see that the ground state energy curve obtained from 
the SAE and CIS approximations overlap with the HF ground 
state. Further, these approximations are unable to reproduce 
the correct dissociation limit due to the restricted HF ansatz 
used in this work, which prohibits dissociation into two open 
shell systems. For the GAS methods, the wave function is 
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written as 

|^GAsW)=CoW|$ 0 )+E C '“WI $ i) 

ia 

+s £<#(*) !*«> + ••• • ( 15 > 

ijab 

Let us consider the CAS* (2, 2) model with two active orbitals 
in the lowest GAS partition. It means, the summation over 
the core orbitals in Eq. (15) runs from one to two. Similarly, 
in the CAS* (2, 3) model we have three active orbitals and so 
on. With the CAS* (2, 2) model Fig. 4 shows that the ground 
state energy curve improves significantly over the SAE and 
CIS approximations and it tends to converge towards the exact 
TDSE ground state. We observe that the CAS*(2,4) model 
with four active orbitals in the lowest GAS partition produces 
a converged result. In order to check the convergence, we fur¬ 
ther increased the number of active orbitals to six and then to 
eleven and found that the CAS* (2,4) model is already con¬ 
verged. This is due to the multi-reference character of the 
CAS wave functions, which allows for the break up into two 
H atoms with singly-occupied orbitals. This is in contrast to 
the underlying restricted HF method, where orbitals are al¬ 
ways doubly occupied. 

After the preparation of the H 2 -like molecule in its ground 
state by imaginary time propagation, we use the 800 nm 
laser pulse of Fig. 3 with tfcEP = 0 and two different field 
strengths, (i) F 0 = 0.05 (8.75xl0 14 W/cm 2 ) and (ii) F 0 = 
0.15 (7.87xl0 15 W/cm 2 ), to ionize the molecule. Since the 
H 2 model is homonuclear, the <pcep = tt field of Fig. 3(b) 
gives the same results for the ionization yield as (ficsp = 0. 
In Fig. 5 we show the ionization yield as a function of the 
inter-nuclear separation. In Fig. 5(a) we compare the results 
between the SAE approximation and TDSE calculations. We 
have scaled down the SAE results by an order of magnitude. 
With the increase of the inter-nuclear distance, the ionization 
yield enhances in both models and then it gradually decreases. 
In the SAE approximation the El peak is at R = 4.1 which is 
well below the exact TDSE result of R = 4.7. An earlier cal¬ 
culation on a similar ID model observed the El peak around 
the same inter-nuclear distance [17]. The SAE approxima¬ 
tion captures the overall behavior of El but it overestimates 
largely the magnitude and peaks at the wrong position and, 
therefore, demonstrates the failure of the model to describe 
the El phenomena quantitatively. The comparison between 
the CIS approximation and the TDSE calculation is shown 
Fig. 5(b). Similarly to SAE, the CIS yield is multiplied by 
0.1. The CIS approximation is unable to capture the feature 
of El with the laser parameters used in the present calcula¬ 
tion. It can describe the features qualitatively at the equilib¬ 
rium separation, but the El phenomena which occurs at large 
inter-nuclear separation can not be described. 

In order to understand the electron correlation effects we 
turn to the discussion of Fig. 5(c). The lower curve corre¬ 
sponds to the CAS* (2, 2) model, i.e., two active orbitals in 
the lowest GAS partition, and one can see that this model pre¬ 
dicts the correct El behavior but the magnitude is much less 
than the corresponding TDSE value. Unlike the SAE and CIS 
approximations, the CAS* (2,2) model nevertheless captures 



R (a.u.) 

FIG. 5. Ionization yield vs inter-nuclear distance R for our model 
H 2 with Fo = 0.05 for (a) SAE approximation, (b) CIS approxima¬ 
tion, and (c) GAS approximations. The CAS* (2, 2) model represents 
two active orbitals in the lowest GAS partition and similarly, e.g, the 
CAS* (2,11) model incorporates eleven active orbitals in the lowest 
GAS partition. The squares are fully correlated TDSE reference data. 


the correct position of the El peak. By successively increasing 
the number of orbitals [CAS* (2,3) to CAS* (2,11) in Fig. 5], 
we achieve convergence of the ionization yield and observe 
perfect agreement with the TDSE calculations (squares). We 
find that at least six active orbitals in the GAS partition are 
needed to describe the El phenomena quantitatively. 

In the next step, we increase the laser intensity to Fq = 0.15 
and observe that the CAS* (2,6) model with six active orbitals 
is again the smallest GAS partition where the results converge 
to the TDSE results, see Fig. 6. The idea of introducing the 
GAS concept is to mitigate the computational constraint aris¬ 
ing due to exponential scaling. By comparing with the exact 
TDSE results, we found, as expected, an enormous gain in 
terms of computational costs. To give an example, the con¬ 
verged CAS* (2,6) calculations took 35 CPU hours on a 2.8 
GHz Intel Ivy bridge CPU while the TDSE needed about 1400 
CPU hours on an AMD Opteron CPU. Albeit the drastic re¬ 
duction in CPU time, the TD-GASCI method with a limited 
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FIG. 6. The same as Fig. 5, but for Fo = 0.15. 


number of active orbitals is still able to describe quantitatively 
the El phenomena, including correlation. 


B. Four-electron LiFl-like molecule 

We now turn our attention to the four-electron heteronu- 
clear LiH molecule for which TDSE simulations of El are in¬ 
feasible. LiH lacks the reflection symmetry of H 2 at the origin 
and thus strong CEP for few-cycle pulses can be expected. 

We choose the geometric center as the origin of the coordi¬ 
nate system and choose Z i = 1 and Z 2 = 3 in Eq. (3), i.e., 
the Li atom is at —R/2 and the H atom is placed at +R/2. We 
prepare the LiH-like molecule in its ground state by imaginary 
time propagation. All numerical parameters are the same as 
for H 2 , see Sec. Ill A. 


1. Ground-state energies 

The ground state energy curve for LiH is shown in Fig. 7. 
In Panel (a) the HF, SAE and CIS approximations with 



2 3 4 5 6 7 


R (a.u.) 

FIG. 7. (color online). Binding energy of the LiH model, (a): 
HF, SAE and CIS approximations together with GAS calculations 
[CAS* (2, •)] with two active electrons, (b): the same but with four 
active electrons in the CAS [CAS*(4, •)]. 


CAS*(2, •) approximations involving two active electrons are 
shown. Similar to H 2 , the multi-reference description is es¬ 
sential and SAE and CIS fail to reproduce the correct disso¬ 
ciation limit. We find that a CAS* (2, 5) is sufficient for two 
active electrons, as seen by comparison with the CAS* (2,10) 
calculation. To test the frozen-core approximation, we use 
four active electrons and obtain the converged result for the 
ground state energy, as shown in panel (b). As expected, 
the energy is lowered but the equilibrium distance is around 
the same internuclear distance. Similar to Ref. [45], the 
CAS*(4, 5) with five active orbitals produces the converged 
ground state energy curve. 


2. Enhanced ionization 

In order to understand the correlation effects in El phenom¬ 
ena, we expose the LiH molecule to a laser pulse with field 
strength Fq = 0.025 (2.18xl0 13 W/cm 2 ) for ^?cep = 0 and 
V’cep = 7T [Fig- 3], This choice is equivalent to the two pos¬ 
sible orientations in ID of the molecule w.r.t. the dominant 
half-cycle. In Figs. 8(a-c), we show the El calculation with 
^cep = 0. We compare the results of the SAE, CIS approxi- 
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FIG. 8. Ionization yield vs inter-nuclear distance R for LiH with Fo = 0.025. In panels (a-c), ip cep = 0, for (d-f) tpcEP = 7r (see Fig. 3). 
(a) and (d) Comparison between SAE, CIS approximations and fully converged GAS methods with two and four active electrons, (b) and 
(e) Comparison between different GAS methods with two active electrons, (c) and (f) Comparison between different GAS methods with four 
active electrons. 


mations and two GAS methods with two and four active elec¬ 
trons in Fig. 8(a). The SAE approximation shows the El peak 
at R = 5.2. Like in the case of H 2 , the SAE approximation 
predicts the El phenomena for LiH, but again the CIS approx¬ 
imation fails to describe the El phenomena. In Fig. 8(b), we 
show different GAS calculations with two active outer elec¬ 
trons. The CAS* (2, 2) model does not predict the correct El 
peak position and it overestimates the magnitude of the ioniza¬ 
tion yield. Further increasing the number of active orbitals to 
three we found that the El peak is shifted towards larger inter- 
nuclear separation. We systematically increased the number 
of active orbitals and found that we need at least five active 
orbitals to obtain converged results for El. The converged El 


peak is at R = 6.1 with the CAS*(2, 5) model. To check 
this convergence we increase the number of active orbitals to 
ten and the observed yield shows a trend of convergence. In 
Fig. 8(c) we show the results with four active electrons in the 
GAS partition and the CAS* (4, 5) model with five active or¬ 
bitals in the GAS partition leads to a converged result. The 
CAS*(4, 5) model predicts the same qualitative behavior like 
the CAS* (2, 5) model. The El peak in the CAS*(2, 5) and the 
CAS*(4, 5) models have the same position but the magnitude 
is different. In terms of computational cost the CAS* (4, 5) 
model is much more expensive than CAS* (2, 5). To give an 
example the CAS*(4, 5) model took 530 CPU hours com¬ 
pared to 24 hours of the CAS* (2, 5) model on a 2.8 GHz Intel 
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FIG. 9. (color online). Sketch of the strong-field ionization scenarios 
of the LiH model. Shown are the binding potential (black solid lines) 
and the modified potential (red dotted lines) by the electrical field 
(green dashed lines) at the peak of the single-cycle pusle F(t) along 
with the valence orbital (black dashed dotted). The orbital energy is 
indicated by the blue horizontal line. Panel (a) is for </?cep = 0, (b) 
for </jcep = 7r (see Fig. 3). 


Ivy bridge CPU. Thus using only two active electrons in the 
GAS partition gives us more than a 20 times reduction in the 
computational cost, a tremendous gain. 

In Figs. 8(d-f), we show the El calculation with lpcep = 7r, 
and the conclusions w.r.t. convergence are the same as for 
Figs. 8 (a-c). As we incorporate more correlation the yield is 
converged with five active orbitals. In Ref. [45] it was shown 
that the ionization yield at equilibrium is larger for tpcsp = 0 
compared to pcep = 7r . We observe that the ionization yield 
with </>cep = 7r is consistently greater than with pcep = 0 at 
the critical inter-nuclear separation in the different GAS cal¬ 
culations. It certainly illustrates the importance of correlation 
contributions in the direction of electron emission from a po¬ 
lar molecule like LiH. 


C. CEP effects in El of LiH 

In this section, we discuss the CEP effects in EL The im¬ 
portance of the CEP of the few-cycle pulse in El of po¬ 
lar molecules was demonstrated on a one-electron system, 
HeH 2+ in Ref. [21, 22] and on a two-electron system, HeH + 
in Ref. [23]. In this work we illustrate the importance of CEP 
effects in El in the four-electron LiH model. As alluded in 
the beginning of Sec. Ill, to observe such CEP effects experi¬ 
mentally would require a pump-probe scheme. The LiH-like 
molecule can have two orientations w.r.t the linearly polar¬ 
ized laser field. In Fig. 9, we show a schematic diagram of 
LiH interacting with a laser field. The interaction with F(t) 
at the peak of the pulse with ip cep = 0 [Fig. 3(a)] is illus¬ 
trated in Fig. 9(a) and with <^cep = 7r [Fig. 3(b)] in Fig. 9(b), 
together with the field-free valence HF orbital. To study the 
CEP effects, the intensities are chosen to illustrate the correla- 



FIG. 10. (color online). Ionization yield vs inter-nuclear distance R 
for the LiH model in the SAE approximation. See Fig. 3. for the 
definition of ip cep- The panel with Fq = 0.005 is multiplied by a 
factor of ten to scale. 


tion effects in the El phenomena in the tunneling as well as in 
the over-barrier-ionization regime. First we use Fq = 0.005, 
0.008 which are below the over-barrier-field strength and rep¬ 
resents the tunneling regime and the others with F 0 = 0.020, 
up to 0.045, which are above the over-barrier-field strengths. 

The CEP effects in the El are shown in Figs. 10 (SAE) and 
11 [correlated CAS* (2, 5)]. In the tunneling regime, the El 
with pcep = 0 shows a completely different trend than with 
the ipcEP = tt- The El peak is well captured with <pcep = 0, 
but with tpcEP = 7r it is absent at the expected inter-nuclear 
separation. The same feature is observed with the CAS* (2, 5) 
model in Fig. 11. The correlated model predicts a larger ion¬ 
ization yield at the critical inter-nuclear separation than the 
SAE approximation. It is important to emphasize that the 
El peak is shifted towards larger inter-nuclear separation with 
Pcep = tt. This is a CEP effect because both, the SAE and the 
CAS* (2, 5) approximations, predicts the same trend. As we 
increase the field strength to Fq = 0.025, the SAE approx¬ 
imation predicts the El peak for both pcep at similar inter- 
nuclear separation. However, in the SAE approximation there 
is a cross-over between the El peak with two different pcep- 
As we increase the intensity, the SAE approximation predicts 
more ionization yield with pcep = tt than with <^cep = 0 
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R (a.u.) 


R (a.u.) 


FIG. 11. (color online). Ionization yield vs inter-nuclear distance 
R for the LiH model in the CAS*(2, 5) model. See Fig. 3. for the 
definition of ipcEP- The panel with Fq = 0.005 is multiplied by a 
factor of ten to scale. 


at the critical inter-nuclear separation. The correlated model 
gives a completely different picture. In this model, the ion¬ 
ization yield is larger for ip cep = 0 than for pcep = tt for 
F 0 > 0.035, which is certainly due to correlation effects. The 
previous calculations with one- and two-electron molecules 
predict larger ionization yield with the pcep = 0 than with 
<Pcep = 7T at the critical inter-nuclear separation [21-23]. 
However, our studies show that there exists a cross-over be¬ 
tween the ip cep = 0 and pcep = tt cases in the over-barrier 
regime as well as the pcep = ~ predicts more ionization yield 
than with pcvv = 0 in the tunneling regime. Therefore the 
CEP of a few-cycle pulse plays an important role in describ¬ 
ing the direction of electron ejection in both tunneling and 
over-barrier-ionization, and most importantly the electron cor¬ 


relation can completely reverse the trend of the SAE approxi¬ 
mation as seen from, e.g, the Fq = 0.035 case. 


IV. SUMMARY AND CONCLUSION 

In this work, we investigated the importance of electron 
correlation effects in enhanced ionization of model diatomic 
molecules by using the TD-GASCI method. The strength of 
the TD-GASCI method is to incorporate the electron corre¬ 
lation systematically. We used ID models of H 2 and LiH 
to study the electron correlation effects. We used the imag¬ 
inary time propagation to calculate the ground state of the 
molecules. We then gave a detailed analysis of El with dif¬ 
ferent laser intensities and illustrated the features that arise be¬ 
cause of electron correlation by using the TD-GASCI method. 
We found that correlation is very important to describe the El 
behavior in the H 2 model. The SAE approximation failed to 
describe El in H 2 model and we would expect the same trends 
in real 3D systems. Another important conclusion from our 
study is the failure of the CIS approximation to correctly de¬ 
scribe the El phenomena. On the other hand the different GAS 
approximations predicted the correct behavior of El and as we 
incorporated more and more correlation in the H 2 model the 
result converged towards the TDSE calculations. For the LiH 
model we found that El persists and that the correlation shifts 
the peak of the El towards larger inter-nuclear separation. We 
also estimated the computational cost to solve the TD-GASCI 
equations with different GAS partitioning and found that it is 
computationally feasible and inexpensive compared to TDSE 
calculations. We emphasized the advantage of two active elec¬ 
trons over four active electrons in the GAS partition for the 
particular problem of strong-held ionization where the va¬ 
lence electrons play a dominant role. Compared to the exact 
four electrons in the GAS partition, the two active electron 
model predicts the yield correctly and with much less compu¬ 
tational cost. The CEP of a few-cycle pulse plays an important 
role in describing the El in the asymmetric LiH model. The 
intriguing feature is the cross-over between the <^cep = 0 and 
Pcep = 7r in the tunneling as well as in over-barrier-ionization 
regime. In conclusion the electron correlation plays a crucial 
role in the El phenomena. 
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